CfPA-94-TH-28 
astro-ph/9406060 

QUOTING EXPERIMENTAL INFORMATION 

MARTIN WHITE & DOUGLAS SCOTT 

Center for Particle Astrophysics, University of California, 
Berkeley, CA 94720-7304 

ABSTRACT 

There was a question raised at this meeting as to the best way to present the 
results of experiments measuring anisotropics in the CMB. Here we will make some 
simple comments about the 3 main competing methods and some suggestions. In 
particular we will give an over-simplified but hopefully useful method of comparing 
GACF numbers with other power spectra. 

Subject headings: cosmic background radiation — cosmology: theories and obser- 
vations — GACF's: justified 



1. Introduction 

There are a number of ways of quoting the level of anisotropy measured by a 
CMB experiment. Some of the most common are: the rms temperature measured by 
the experiment T rms ; the amplitude and correlation angle of a Gaussian Autocorre- 
lation Function (GACF) Co and 9 C ; the amplitude of a 'flat' or Harrison-Zel'dovich 
power spectrum; and the amplitude of a 'standard' spectrum such as CDM. 

The rms temperature measured by the experiment is the quantity most directly 
related to the data, but also the least informative when it comes to comparing 
different experiments, since it depends on the window function and calibration pro- 
cedures etc. of the experiments. It also does not include the sample 1 and cosmic 
variance associated with the measurement. All the other methods choose to quote 
the amplitude of fluctuations by assuming the power spectrum has a certain form 
and fitting for the amplitude. This quantity is independent of the window func- 
tion normalization, calibration etc., and if the amplitude really is fit properly to 
the data, cosmic and sample variance are included in the result automatically. As 
long as the same form is assumed for all experiments, such numbers can be safely 
compared from experiment to experiment and with theories. 

Since the ability to compare experiments is certainly something we would all 
like to have, it seems then that fitting a power spectrum is the way to go. One 
suggestion (Scott Dodelson, this meeting), is to define a 'standard' model such as 
CDM with a fixed Q# and h, and to quote the amplitude of such a spectrum. This 
has the advantage that, to the extent that such a model correctly describes the 
fluctuations on the sky, the amplitudes so obtained can be trivially compared from 



experiment to experiment: they should be all the same. The drawbacks of this 
particular form for the power spectrum are that it involves an enormous amount of 
theoretical prejudice at the level even of quoting the experimental measurements, 
and that the 'standard' power spectrum has a large amount of structure in it. 
Should the model turn out to be something other than 'standard CDM' or if or 
h turn out to be different than was assumed, this structure has to be deconvolved 
from the measurements before a new power spectrum is used. We note in passing 
that comparison of C/s from different codes has not been done in detail to let us 
understand the reproducability of the curves. 

At the other extreme are the flat spectrum and the GACF. These are not at all 
well motivated as models of our sky, but they are extremely simple power spectra. 
While measurements on different scales are not expected to give the same ampli- 
tude (Doppler peaks give a larger amplitude for experiments on degree scales, for 
example), it is not difficult to remember a few numbers at a range of scales for 
comparison, and they have even been tabulated for CDM 2 . 

So what are the merits and disadvantages of GACF's and flat spectra? Histor- 
ically experimental results (upper limits until recently) have been quoted in terms 
of the amplitude of a 'most sensitive' GACF. This method pre-dates the emergence 
of window functions into the popular consciousness and gives an alternate view of 
the 'scale' at which the experiment is sensitive and the degree of the sensitivity 3 ' 4 . 
As discussed in references 2 and 5 the GACF analysis is straightforward to un- 
derstand in terms of window functions, and a conversion from GACF numbers to 
'flat spectrum' numbers is easy to do 2 ' 6 ' 7 ' 8 . Note that the advantage of having a 
readily understandable scale in your fit (9 C ) is also a disadvantage when it comes 
to comparing experiments, since values of Cq for different (arbitrary) 9 C cannot be 
directly compared. The flat spectrum has no such problem since it is featureless 
and easily understandable, but does not directly contain the information on the 
window function that the GACF encodes. 

Here we would like to make some basic remarks about the relation of the GACF 
to the window function and about converting from GACF's to flat spectra and 
back. The details of such a conversion have been treated before, but we would like 
to present a simple (and not necessarily very accurate!) method which hopefully 
illustrates the general idea. 

2. To and from the GACF 

Since our aim here is simplicity rather than rigour, we are going to assume that 
all CMB experiments on small scales are two-beam square-wave chop experiments, 
with peak-to-peak chop a and (gaussian) beam width a. Under this approximation 
the window function is of the well known form 

W e = 2[l-P e (cosa)}e-^ i+1 ^ 2 

~ \{iafe~^\ 

where in the second line we have made the approximation a < 1, or more cor- 



rectly that £a < 1 for the £ range of interest. This is not a tremendously good 
approximation, but brings out the essential features of what is going on. 

Now in the language of power spectra and window functions, the most important 
quantity measured by an experiment is the power through the window, which is 
given by 

Power = — ]T(2£ + l)C e W e . (2) 

As experiments become more sensitive, other quantities (e.g. correlations) are going 
to play an important role 9 , but for now the power is what most experiments are 
sensitive to 5 . For experiments on small or intermediate angular scales we can ap- 
proximate the sum as an integral and replace (2£ + 1) — > 21. If the window function 
drops off at small £ we are safe in extending the lower limit of the integration to 
£ = 0. The quantity normally plotted with the window function is £(£ + l)Ce £ 2 Ce 
in terms of which 

Power [°° e 2 C e Wt d\og£, (3) 

4-7T J-oo 

which explains the logarithmic £-axis normally used in such plots. 

What do the two 'competing' power spectra look like? The flat spectrum is 
extremely simple 

£(£ + l)C e = —Qlzt = constant. (4) 
5 

The factor of 24-7r/5 is there simply to soak up the constant in Eq. (2) so that Q\ &t 
is the coefficient of W 2 a la COBE. The GACF power spectrum on the other hand 
is 

£(£ + l)C £ = 2nC £ 2 e 2 c e- 1 / 2i2e *. (5) 

Notice now the similarity between £(£ + l)Ct for the GACF and the expression for 
the window function of Eq. (1). Keeping in mind that the power measured by the 
experiment is the fixed quantity defined by the data, it will come as no surprise 
that the 9 C which gives the lowest C is 9 C = \[2o ~ 0.6xFWHM. This is the value 
of 9 C which maximizes the power through the window function for a GACF [using 
Eq. (1) and Eq. (5) in Eq. (3)]. All that varying 9 C in the fit is doing is matching the 
power spectrum to the window function of the experiment, and this is another way 
of encoding the information about which ts the experiment measures. In the highly 
simplified approximation for window functions given by Eq. (1), the correlation angle 
of the experiment is set by the beam size, which also sets the peak of the window 
function, which is £q = 1/a. In practice this last relation is very poor (see Table 1), 
since getting the peak correct requires accurately modelling the left hand rise of 
We well for 'large' £a unless, a 3> a. This needs a detailed consideration of each 
experiment 10 . 

We can do some easy integrals to see how the GACF-to-flat spectrum conversion 
works in our over-simplified but hopefully instructive approximation. The power 
through the window of the flat spectrum is proportional to the area under the 



window function vs. log £, or 

6 f°° 

Power Q| / di 2 r 2 W £ 
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While the power through the window for a GACF is 



Power = -C 2 c / 
2 Jo 
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where in the second line we have evaluated the integral for 6 C = \plo. Relating these 

1 12 

two we see that Qaat — 0.46Cg , which compares well with numbers obtained using 
the accurate Wis and the 9 C ' 's quoted for each experiment (see Table 1). 
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Table 1: Summary of the peak of the window function and the C to Qfl a t conversion factor 
for current CMB experiments on small and intermediate scales. Note that the peak of the window 
function does not correspond very well to l/er, showing that <j>ctis not a good approximation for 
most experiments. However, the C]J 2 to <3 flat conversion is close to 0.46-0.50, for many experiments, 
as predicted by simple calculation. 



Experiment 


to 


1/(7 


e c 


Q/c 1 / 2 


Tenerife 


20 


25 


4°.0 


0.50 


SP91 


66 


96 


1°.5 


0.44 


SK93 


71 


93 


1°.2 


0.49 


Python 


73 


180 


1°.0 


0.47 


ARGO 


107 


156 


0°.5 


0.42 


MAX 


158 


270 


0°.5 


0.45 


MSAM2 


143 


289 


0°.5 


0.44 


MSAM3 


249 


289 


0°.3 


0.50 



We can do a little better for some experiments by noticing that they perform 
a double-difference ('triple-beam') rather than a single difference ('double-beam') 
measurement. This means that the window function is better represented by the 
functional form Wg oc g 4 e - £ ( i + 1 ) a ~ than by the lower order expression in Eq. (1). 
This changes the peak to £q = \[2ja^ and the power of a GACF through this window 
peaks at 9 C = a. The comparison of a GACF and a flat spectrum then leads to 
Qflat/C 1/2 = 2^/9 ~ 0.50. We would expect Tenerife, SK93, Python and MSAM3 
to be better characterized by this double-difference approximation. A glance at 
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Table 1 shows that the Q/C ' ratios are indeed higher. 



3. Conclusions 



In reality the approximations made above are not quite adequate for accurately 



converting from GACF numbers to other measures of the power, but the procedure 
has the virtue of being very straightforward: one calculates the power in Eq. (2) for 
whatever power spectrum was used in the quoted answer. Assuming that this is the 
important information that the experiment provides, you adjust the normalization 
of your favourite power spectrum to match this number. Note that as long as you 
use the same window function in both calculations, the Wg normalization doesn't 
matter. This is a boon because getting the shape of the window function is usually 
a lot easier than determining the height! Scaling from numbers quoted for power 
spectrum fits has the advantage that all of the experimental and theoretical errors 
will be included in such a fit. 

Hopefully working through this over-simplified example has given some insight 
into the relation between the well loved GACF method and the £-space measures 
preferred by theorists. 

The bottom line: If you want a quick and dirty Co to Q conversion, then a factor of 
two is good at the 10% level for all experiments to date. 
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